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Theory of Ultrasonic Dispersion in Local Phonon Systems Coupled with Conduction 
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The physical origin of frequency dependence in elastic constants, which are often found in an 
ultrasound propagation in filled skutterudites and clathrate compounds, is investigated theo- 
retically. This dependence arises from a coupling between the acoustic phonon and some optical 
phonons, which strongly interact with electrons. Using a self-consistent ladder approximation 
together with a pseudofermion mapping of the phonon to the single site Holstein Anderson 
model, a soft mode of the optical phonon at zero frequency is shown to emerge. The temperature 
dependence of the spectral weight of this soft mode shows an activation-type behavior, which 
is characterized by the optical phonon frequency. These features can generate the frequency 
dependence and the shoulder in the elastic constants observed in some filled skutterudites and 
clathrate compounds. 

KEYWORDS: skutterudite, clathrate, ultrasonic dispersion, elastic constant, local phonon, Holstein An- 
derson model 



1. Introduction 

Recently, the frequency dependence of ultrasonic ve- 
locity in some filled skutterudites and clathrate com- 
pounds has attracted a great deal of attention. In 
these systems, it is believed that a rare-earth ion in a 
rigid cage moves rather independently in the anharmonic 
potential of off-center minima. In PrOs4Sbi2, for exam- 
ple, it is expected that these ionic degrees of freedom 
are important to heavy fermion superconductivity and 
its large effective quasiparticle mass.^ The data of elas- 
tic constants (EC) in these compounds exhibit frequency 
dependence at temperatures ranging from 10 to 40 K. 
This temperature range is on the order of the Einstein- 
like modes of the rare-earth oscillation in these com- 
pounds.''"^^ For PrOs4Sbi2, large softening of acoustic 
phonons is also reported. Recently, Iwasa et al., ob- 
served the development of a quasi-elastic peak at low 
temperatures in Pr-based filled skutterudites.^"^ 

Usually, the frequency dependence of EC manifests as 
a shoulder in its temperature dependence. This is well 
fitted by a Debye-type formula: 

C(z) = C.o + ^^y%, (1) 

1 + (ZT)^ 

where C(z) is the frequency-dependent EC and C'oo = 
C(oo) and Co = C(0) are constants in frequencies. In 
the present paper, we use z as the frequency of dy- 
namic quantities. The lifetime t in eq. (1) is given as 
T = To exp(i?A/r), where Ea is a characteristic en- 
ergy scale and T is the temperature. The exponential 
T-dependence gives the shoulder at zr ~ 1, at approx- 
imately r ~ 40 K for Ea 200 K, even though the 
frequency of the ultrasound is of the order of mK. This 
is called "ultrasonic dispersion" (UD). This form of the 
relaxation time appears to be related to a kind of ther- 
mal activation-type process, so that in ultrasonic exper- 
iments it has been interpreted that there is an off-center 



potential and the origin of the relaxation is due to the 
thermal hopping between the off-center sites. In this re- 
spect, we discussed a possible scenario for the realization 
of heavy fermions in SmOs4Sbi2^^ in the strong coupling 
limit of electron-phonon couplings, ' which is a nat- 
ural extension of a two- level system. ^^"^^ Since the UD 
is observed only in a specific mode, it is claimed that 
the off-center degrees of freedom have a degenerate and 
anisotropic ground state. On the other hand, the results 
of a neutron scattering experiment in PrOs4Sbi2 suggests 
that there are no anisotropic charge distributions of the 
Pr-nuclei at low temperatures in particular. 

In order to discuss the frequency dependence in EC, 
it is convenient to treat the two frequency regions sepa- 
rately. One typical region is zr ^ 1, where the time scale 
of the ultrasound is much longer than that of the relax- 
ation mode. At this stage, the origin of the relaxation 
mode is unknown. The relaxation mode is scattered in 
a much shorter time than the oscillation period of the 
ultrasound. Thus, the strain caused by the ultrasound 
can be regarded as static. The ultrasound attenuation 
due to the anharmonicity of the lattice in this region 
was discussed by Akhieser,^"^ and later by Woodruff and 
Ehrenreich,^'' using the semi-classical approach (Boltz- 
mann equation). For the case of zr ^ 1, the Landau and 
Rumer theory, which is "golden rule" treatment, can be 
applied. Although the former description appears to be 
good even in zr > 1 qualitatively, there is no theoretical 
reasoning to use the former in this region. 

Using the Woodruff and Ehrenreich theory, the fre- 
quency dependence in the EC was discussed in KTaOa by 
Barrett^^ from both theoretical and experimental points 
of view. In KTaOa, the activation factor exp(a;o/2^) also 
appears in the ultrasound velocity, where loq roughly 
corresponds to the frequency of the optical mode. This 
factor arises from the specific heat of the free opti- 
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cal mode through the time-dependent Boltzmann equa- 
tion. In comparison, the theory for the frequency depen- 
dence in filled skutterudite and clathrate compounds is 
more difficult, since we need to bridge the above two re- 
gions continuously and to take into account the electron- 
phonon couplings. 

In the present paper, we discuss a simple interpreta- 
tion of this problem. In §2, we review the Green function 
formulation in electron-phonon systems. We then discuss 
two types of couplings between the optical and acoustic 
phonons, and give the phenomenological interpretation 
of the Debye formula used in the ultrasonic experiments. 
In §3. wc investigate the Holstein- Anderson model at 
finite temperatures using the self-consistent ladder ap- 
proximation. We then show the results of the renormal- 
ized sound velocity, giving the dispersion in EC, together 
with the nature of soft optical phonon modes. In §4, we 
discuss the applicability of the present theory to the real 
materials and future problems. Finally, in §5, we sum- 
marize the results. 

2. Coupling between acoustic and optical 
phonons 

In this section, we discuss how eq. (1) is derived in 
interacting electron-phonon systems. Note that the dis- 
placements due to the acoustic modes near q = include 
contributions from optical modes at q = 0. This essential 
point is described in §2.1. In §2.2, we discuss the case in 
which the anharmonicity of the lattice is important. 

2.1 Harmonic coupling 

The classical Hamiltonian for the phonon system in the 
harmonic approximation is written using a force constant 
matrix K as 



ph 



E E ^2 



o i,j,n,m in jmi 



(2) 

where P^iX^) is the momentum (displacement) variable 
of the ion in the ji direction located at the n-th site in 
the i-th unit cell. M„ is the mass of the ion at the n-th 
site in a unit cell. Introducing variables = Pl^ /\/M^, 

^ VM^Xt^ and k^^l^^^ ^ VM^'' 
the Hamiltonian (2) is reduced to 



= E^^rnPL+E E E<I 

\i,n,i^ iJ tJ,iy=x,y,z n,m 

= ^EfE^^"(q)K(-q) 

q \ n,iy 



\,n,m m jm 



+ 



E E^"rm(q)<(q)^™(-q) 



(3) 



fj.iy—x.y,z n,m 



where we have introduced the Fourier components of 
each variable in the second line. After diagonalizing the 
/c-matrix, we obtain all of the eigenvalues in principle as 
in textbooks.^* In the usual step of moving to the corre- 
sponding quantum Hamiltonian, we use the diagonalized 



basis. Although we can diagonalize the /c-matrix and ob- 
tain the energy eigenvalues of every phonon branch, we 
here introduce even and odd parts, and k°, with re- 
spect to q, and concentrate on the low q limit of the 
phonon system. In order to determine the phonon veloc- 
ity, we need only retain up to the second-order contribu- 
tion in q. Thus, we write k'' = fc^"-* -I- fc*^^^ and k° = k^^\ 
Dividing k into k^ and k° is meaningful only in the region 
near q = 0. Near q = and z = {z being the frequency 
of phonons), k° = 0(|q|) acts as the hybridization be- 
tween some optical and acoustic modes at q ~ 0. Since 
we are interested in the sound velocity of the system un- 
der which optical phonons strongly interact with conduc- 
tion electrons (this comes from the fact that UD has an 
activation type relaxation) , dividing k into the even and 
odd parts is meaningful. For large |q| > {a being the 
lattice constant), this characterization loses its meaning, 
and k° should be regarded simply as the hybridization. 

In order to discuss the quantum mechanical Hamilto- 
nian, we first diagonalize fc*. We define the creation and 
annihilation operators of the mode A with the energy 
eigenvalue ojqx (related to the eigenvalue of k^, not k) as 
follows: 



1 



(4) 



Here, Xqx is the displacement operator that diagonalizes 
fe*^. Using this set of displacement variables, Hamiltonian 
(3) can be written as 

Hph = E'^I^^^A^q^ + J) + J E ^qAfc^A'(q)*-qA'- (5) 



qA 



qAA' 



Next, we introduce the retarded Green function Dqxx'{t) 
in the usual manner: 

^qAA'(i) = -l([iqA(t),i-qA'(0)])^?(t), (6) 

= ^Jdze-^^'DqXX'iz). (7) 

The equation of motion for D^xx'it) has the following 
form: 



E 



-'^qA)^AA" - 



fcA"A(q) + ^AA"(-q) 



xD^x"X'it)=S{t)6xx'- (8) 
It is useful to employ Dqxx' (t) , which is defined as 

(9) 



DqXX'it) ^ V2^^qAA'(i)V2^ 



qA' 



The Fourier components of Dqxx'{t) obey the following 
coupled equations: 



E 



'qA 



(^^ - ^qx)^)^X" - «AA"(q) 



1 



-v/2WqA" 

xDqX"X'{z)=6xx'- (10) 

Here, we have written the mixing term ^{ky,^{q) + 
^AA"(^q)) '*AA"(q)- The energy eigenvalues are deter- 
mined by the poles of Dqxx{z)- It is evident that these 
poles coincide with those in usual basis, which diagonal- 
ize fc« + k°. 
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If we take the electron-phonon coupling into account, 
we must introduce the polarization functions IlqAA" (z) 
into the left-hand side of eq. (10). Defining {Dq^{z))x\/> 
as the matrix in the large brackets of the left-hand side 
in eq. (10) divided by ^y2uJqX\/^^^q>J'^ the equation of 
motion for Dq\iix> (z) in the interacting electron-phonon 
system can be written as 



E 



In order to clarify the importance of kxx" (q) , it is in- 
structive to consider the simplest case in which there 
exist one acoustic (its energy (Dqi) and one optical mode 
('i'q2) in a system we consider. Furthermore, the electron- 
phonon couplings are nonzero only for the optical mode 
and KAA(q) = by definition. Here, we use the words 
"acoustic" and "optical" with respect to those in the ba- 
sis Xqx- In the diagonalized basis, there is no mixing be- 
tween the acoustic and optical phonons in the harmonic 
Hamiltonian (2). Then, eq. (11) is reduced to 



-Ki2(q)i 




- l^q2 - 2a;q2nq22(z) 



-Dqll(z) Dql2{z) \ _ ( 2cZ'ql 



-Dq2l(2:) Dq22{z) 



(12) 



2c:;q2 ) 

We find the form of the Green function for the acoustic 
mode as 

''^l2(q)P 



qi 



/"•.2 
^q2 



2u)q2^q22{zy 



(13) 



Noted that in eq. (13), Ki2(q) — Ko(9)|q|j <^q2 — ^ — 
const, and cDqi ~ Mo(9)|q| for small q, where q = q/|q|. 
Although the expression a)qi ~ uo(g)|q| is not satisfied in 
general, it is a reasonable assumption for filled skutteru- 
dite and clathrate compounds. Diagrammatically, eq. 
(13) is represented by Fig. 1(a). Thus, the frequency- 
dependent phonon velocity m(z, q) is given formally by 

l«o(9)P 



u {z,q) = u^iq) + 



z2-02-2f]n^22(^)' 



(14) 



In the limit of ultrasound, we can ignore z <^ ^ and set 



n 



ij22 



— iVqZ in eq. (14), leading to 



v?{zA) 



l«o(9)P 



2VlIi'^ 1 



i2Tqz/{^- 



u^{Qo,q) 



u^{0,q) -u^{(X),q) 



1 



2U'.y 



(15) 



appears around zero frequency, is related to the many- 
body effects between the optical phonon and conduction 
electrons). In the above simple model, we have ignored 
Ilqnm with (n, m) ^ (2,2). This is because we consider 
that most of the important properties come from the 
self-energy of the optical phonon. If these come from the 
degrees of freedom in the lattice vibrations, it is clear 
that i?A cannot appear without the finite energy excita- 
tion of the optical phonon. 

2.2 Anharmonic coupling 

The discussion in the previous subsection, however, is 
not directly applicable to the rare-earth mode in filled 
skutterudites because the optical phonons that hybridize 
with the acoustic phonons near q = are those with 
even parity in that system. The rare-earth mode cou- 
ples with the acoustic phonon only through the anhar- 
monic coupling so long as we ignore the electron-acoustic 
phonon interactions, as shown in Fig. 1(b). Since the 
ultrasound generates strain £r(q) with the wave num- 
ber q and the symmetry F = F^, FJ and F^ in the 
cubic Oh point group, the rare-earth mode can couple 
with these strain-fields by constructing the direct prod- 



uct: F4 (g) F4 = F^ © FJ © F^. ((rr(q) is expressed by 



the acoustic phonon fields.) Thus, we have coupling due 
to the anharmonicity as 



E-9r £r{q){(t)-(l)-)-qr, 



(16) 



where 4>- is the phonon field with the odd parity and 
is the coupling constant. Here we use the symbolic 
notation ((/)_<^_)_qr. For instance, {(j)-4>-)_ should 

be read as the Fourier component of {(pixfj^iy) with the 
wave number — q. Corresponding to ^ the coupling 
between the optical and acoustic phonons in the previous 
subsection can be rewritten as 



V+ = ^5rit£r(q)5-qr- 



(17) 



The EC's are renormalized by the susceptibility of 
(0_0_)r (denoted as ^(^^{z.fi}). A diagrammatic 
expression is shown in Fig. 1(b). Information of is 



required in order to obtain an explicit form of the sound 
velocity. This is more difficult than in the case of the 
previous subsection. In phenomenological treatment, we 
can assume c< 1/(1 — izTq), which is the same form 
as that in eq. (15). 

We summarize this section as follows: 



where Tg = F^/(f2 -I- 2H'-). Equation (15) corresponds to 
the phenomenological expression eq. (1) used in the anal- 
ysis of the ultrasound experiments. In phenomenologi- 
cal theory, the temperature dependence of Tq is assumed 
to be proportional to exp{E/\/T). The activation scale 
Ea is expected to be of the order of il. The fact that 
the ultrasound experimental results are well explained 
by eq. (1) indicates the existence of a soft mode of the 
phonon (in this simple example, the soft mode, which 



■AA/W\A- 



(a) 




(b) 



Fig. 1. Diagram of the hybridization process between the acous- 
tic and optical phonons near q = 0. The straight (wavy) hne in- 
dicates the acoustic (optical) phonon propagator, (a) Harmonic 
term, (b) Anharmonic term. 
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• High-temperature region: the sound velocity is given 
by "original" sound velocity Uo{q) subtracted by 
the contribution from the "non-interacting" opti- 
cal phonon \Ko{q)\'^ / {Q^ + 2nn'.) (see eq. (15)) be- 
cause of the smallness of ZTq. This corresponds to 
the sound velocity obtained by diagonahzing k"^ + k° 
with renormalized (l. Thus, the sound velocity is ob- 
served as a smaller value than that at lower temper- 
atures as explained below. 

• Low-temperature region: the relaxation time be- 
comes very large around z = 0, and eventually zTq 
exceeds 1. Note that the relaxation time of the opti- 
cal phonon at the original pole is much smaller than 
Tq. As a result, the second term in eq. (15) docs not 
affect the sound velocity. Thus, we observe only the 
"original" sound velocity Uo{q). 

These aspects are similar to those discussed by Yamada 
for the critical fluctuations of lattice systems. In these 
studies, the system is assumed to be located near the 
structural phase transition and the imaginary part of the 
polarization was phenomenologically introduced. In this 
paper, however, we will calculate the explicit tempera- 
ture dependence of Tq, namely oc cxj>{E^/T), based on 
a relevant microscopic model. In the next section, we in- 
vestigate a simple model with respect to the discussions 
in §2.1 for simplicity. However, in principle, it is straight- 
forward to discussion the same topic based on §2.2. 

3. Model and Calculation 

In this section, we concentrate on a model in which 
there is only one electron band and one local Einstein 
phonon. To carry out the complete calculation, it is de- 
sirable to evaluate the properties of the model with two 
kinds of phonoiis (whose Green functions are formally 
given by eq. (12)) and conduction electrons. However, in 
order to obtain a qualitative understanding, we restrict 
ourselves in one local phonon coupled with conduction 
electrons and approximate the frequency dependence of 
the sound velocity in §3.3. 

One of the simplest models is the Holstein model,^^ 
which is given as 

aij i 

+ T.9{4.Ci.-l)ibl + h), (18) 

where represents the hopping of conduction electrons, 
is the creation operator of electrons at site i and spin 
a, 6; is the phonon creation operator at site i with Ein- 
stein energy Qe, and g is the electron-phonon coupling 
constant. This model has been discussed for the past five 
decades. Recently, the phase diagram of this model 
at T = was discussed based on the dynamical mean 
field theory (DMFT)'''' rising the numerical renormaliza- 
tion group (NRG) as an impurity solver. '^^"'^^ The phase 
diagram of a somewhat different model^^ was discussed 
by DMFT, but with exact diagonalization. For finite tem- 
peratures, data is lacking for frequencies smaller than 
the temperature in the NRG method. Recent develop- 
ments in calculating spectral functions in NRG enable us 
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to obtain the spectral functions roughly in the range of 
z >T X 1/10.'^^'''° However, it is insufficient to discuss 
the soft mode at finite temperatures. Thus, in the fol- 
lowing subsections, we introduce a simple self-consistent 
treatment of the phonon system, which can capture the 
essential aspect of the emergence of the soft mode at low 
temperatures. For the connection with the single site ap- 
proximation, which is widely used for experimental anal- 
ysis of the quadrupolar degrees of freedom in f-electron 
systems, '^"'^ we restrict ourselves to the single-site problem 
of a local phonon interacting with conduction electrons. 
This model is called the Holstein- Anderson model, which 
will be explained below. 

3.1 Self- consistent theory in pseudo-fermion represen- 
tation of phonon 
First, we explain a pseudo-fermion mapping of phonon 
operators. The displacement operator X — b'' + b has the 
matrix elements in the phonon Hilbert space as follows: 

Xnm = {n\X\m) 

= \/rn5m,n+i + \/nSm,n-i, (19) 

where n and m are eigenvalues of the number of the 
phonons (6^6). We introduce pseudo-fermions a|j, which 
create the state with phonon numbers n = 0,l,2,---, as 
reported by Abrikosov for the case of "spin" in the Kondo 
problem."'^ Using the pseudo-fermions, the displacement 
operator is represented by 

X = J2aiXnmam. (20) 

In this representation, the Holstein- Anderson model 
with the Coulomb repulsion of the electrons U = be- 
comes 

H = ^[ekcj^^^CkcT + (wckcrcj^ -I- h.c.)] -|- ^ nO_Eat a„ 

kc n 

+ 9 X^^^OaCOcT - ^)alXnmam + ^ A(aJ,a„ - 1|?1) 

(T n 

where cj^(at ) is the on-site conduction electron (pseudo- 
fermion) creation operator with spin ct, ek and v char- 
acterize the conduction electron dispersion and the hy- 
bridization, respectively, and cj^^ is the creation operator 
of the conduction electron with the wave vector k and 
spin a. We have introduced a Legcndrc multiplier A in or- 
der to prohibit double occupancy of the pseudo-fermions. 

Next, we explain a self-consistent treatment of this 
impurity model. The method explained here is simi- 
lar to the self-consistent ladder approximation in the 
Coqblin-Schrieffer model with crystalline-clectric-field 
states of Ce impurities^^''*^ and non-crossing approxima- 
tion (NCA) of Anderson model. ''3. 44 \Ve define the Mat- 
subara Green functions of a„ and cqo- in the imaginary 
time r as follows: 

Anmir) = -{Tran{T)al{0)), (22) 

G{t) = -(T,co.(r)cL(0)), (23) 

where Tt is the time-order operator, and we omit the 
spin dependence in G{t) hereafter. 
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The non- interacting Green functions A"„(ici;„) and 
G'^{iujn), are given as 

where 
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(24) 
(25) 



E 



(26) 



where ujn is the fermionic Matsubara frequency a;„ = 
{2n + l)nT, In order to take into account the electron- 
phonon couphng, we consider the self-energy diagrams 
T,^{iujn) and ^^(ia;) shown in Figs. 2(a) and 2(b), re- 
spectively, as: 



sr(*wn) 



2g2r^G(ze™) 



E 

n'm' 



SG(it^n) 



XgA (*£m — iljJn)XnmXn'm' ■ (28) 



^^^E E 

iem nmn'rn' 



Here, XGA is the solution for the diagrammatic equation 
given in Fig. 2(c), which is reduced to the summation of 
the infinite series of the ladder diagram Xga constituted 
of G and A. This is analytically given by 



E 



[Snn' -sE^Ga"" {i'^l)Xm'n']XGA{m) 



n nm / 

= Xga 



\^n). (29) 

where vi = 2'kIT is a bosonic Matsubara frequency. In 
matrix form, eq. (29) is expressed as 

XGA^[r-gx%AXr'xGA- (30) 
The explicit form of Xga given by 

X^jC'^im) = -Ty^G{ier,i)Anm{m + iern)- (31) 

Note that the first-order term in g in Ha and Eg van- 
ishes, together with the term ^ in the interaction term 
in eq. (21). 

Using E^™ and Eg together with Dyson's equation: 

+ ^ AL,(^c^„)E:^''"'(^c^„)A„,,„(^c^„)(32) 

G(zc^„) = G°(^c^„) + G°(^c^„)EG(^c^„)G(^c^„),(33) 

we can obtain the self-consistent Green functions. In ac- 
tual calculations, it is better to treat the retarded Green 
functions E^™(z -t- iO) and Eg(z -|- iO). Carrying out the 
analytic continuations, we can obtain the set of equa- 
tions for Y7^{z + iO) and Eg (2 + iO), as shown in the 
Appendix (we hereinafter write z -|- iO as simply z). 



mnm'n' 


m'n 
■^GA 




> 




— > 



gXnm gXn'n 

(a) 





m n 



— < — 




gXm'n' 
^>-- 


y n'n 
^GA 


m m 





(c) 



Fig. 2. Self-energy diagrams (a) T,Q(iujn) and (b) S™(jaj„). (c) 
Diagrammatic equation for xga - The solid (dashed) line repre- 
sents G{A). 



The phonon Green function D{t), which is the most 
crucial quantity in the present paper, is given by 

D{t) = -{TrX{r)X{0)) 



^ ^ -^Jim-^n' m' -^mn' ('^)-^m'n ( 



(34) 



where we ignore the vertex corrections in the last line. In 
this approximation, the phonon Green function in Mat- 
subara frequency vi is given by 



D{m) = '^XnmXn'r,l'T'^A„^n^{iuJl, +Wi)Ara'n{i 



nn'mm' 



(35) 

The main feature of the approximation scheme ex- 
plained in this subsection is to take into account "phonon 
fluctuations" . As such, the results obtained in this 
method might overestimate the fluctuations because we 
ignore the vertex corrections in eq. (34). Despite this, the 
present approximation can visualize the essential points 
of the low-energy phenomena in the present electron- 
phonon system. 

3.2 Numerical results 

In this subsection, we show the numerical results of 
the Holstein- Anderson model (21). For the conduction 
electron, we assume the Gaussian density of states p{z): 



1 



^exp(-7^), 



(36) 



which corresponds to an infinite dimensional hyper-cubic 
lattice with the nearest neighbor hopping t. We set 
t — 0.2 and VIe = 0.01 throughout this subsection. The 
hybridization width ImA(z) is given by using p[z) as 

-lmA(z) = 7rw2p(z), (37) 

where w is a constant. The real part of A(z) is calculated 
by the Kramers-Kronig relation 

dx ImA(x) 



ReA(z) = P y" ^ 



(38) 
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where the integral is taken as its principle value. 
In the numerical calculations below, we introduce 
the cutoff number for the pseudo-fermion: TVcut, *-e., 
ao, ai, ••• , aN^,^^^i. Due to this restriction, we can 
calculate the dynamical quantities only in T <, (A'cut — 
1)Qe- However, as shown below, this restriction does not 
become serious in the temperature regions of interest. 



exp(-0.84nE/T) X 0.0074^.-®' 




0.0 0.2 0.4 0.6 0.8 1.0 

Fig. 4. Temperature dependence of ImS(3(0) from Fig. 3. The 
dashed line is the curve of 0.0074 X exp(-0.84S7£;/T). 



- Electron self-energy - 

We show the result of the imaginary part of the elec- 
tron self-energy ImSg in Fig. 3. As the temperatures 
decrease, a dip structure is developed around \z\ < ^Ie- 
This simply represents the fact that an electron with 
an energy below fl^ cannot emit the Einstein phonon 
with energy ^Ie- Thus, the temperature dependence of 
ItoSqIz) for z ~ becomes the activation- type, oc 
exp(Q:/r), where a is a constant on the order of ^Ie, 




Fig. 6. Low-frequency phonon Green function with same param- 
eters as in Fig. 5. 

as shown in Fig. 4. 

- Phonon Green functions - 

In Figs. 5 and 6, we show the frequency and tempera- 
ture dependence of the phonon Green function —lmD{z) 
in high- (Fig. 5) and low- (Fig. 6) frequency regions. The 
position of the main peak around fl e decreases as tem- 
peratures decrease, as shown in Fig. 7. The peak width 
also decreases at low temperatures, making the peak 
sharper. This main peak corresponds to the excitations 
oo <-> ai, ai ^ a2, ■ ■ ■ in the pseudo-fermion picture. 
At low temperatures, below T ^ O.^VIe, a low-lying ex- 
citation becomes prominent, as shown in Fig. 6. This 
low-energy soft mode arises from the development in the 
off-diagonal elements of pseudo-fermion Green functions 
Aoi(O), Ai2(0), • • • , and A„„(0) for n ^ 0. 

Note that the slope at the origin maintains its lin- 
ear frequency dependence, i.e., — ImD(z)/z = const, 
near z = 0, as shown in Fig. 8. Thus, this peak is not 
generated by fictitious increases of the dynamical sus- 
ceptibility lmx{z)/z that frequently arise in NCA at 
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low temperatures.^^''** The temperature dependence of 
this slope is shown in Fig. 9. At intermediate tempera- 
tures (0.1 < T/VIe ^ 1), the temperature dependence of 
the slope is well fitted by the activation-type function, 
exp(a'/r)/T, where a' is a constant on the order of^ls- 
At low temperatures [T/VLe ^ 0.1), the data deviates 
from the exponential dependence and exhibits 1/T^ be- 
havior. These behaviors are discussed below. 

In order to clarify the origin of the temperature de- 
pendence of ImJ5(z), it is useful to consider ImS^™(0). 
In the lowest order ladder diagram of xga, ImE^™ is 
written as 



25^ 



—nB{y)lmxGiy) 

IT 



iy)Xr, 



(39) 



N 

Q 



10' 



10' 



T N 




10~ 



Fig. 9. Temperature dependence of — lim^^o Ii^i'C^)/^ fof vari- 
ous g. The lines are curves of oc cxp(a'/T)/T. The other param- 
eters are as in Fig. 5 



Here, is the Bose function and Imxc is given by 



ImxG(^) 



-[npiy + z) -npiy)] 



X lmG{y + z)lmG{y). 



(40) 



At low temperatures {T < O.IHe in Fig. 9, referred 
to herein as region I), taking into account the fact 
that the temperature dependence in lmiz^o^™-D(z)/ z 
comes from the contribution of the low-energy peak in 
ImAnrm we can replace ImAn'm' in eq. (39) by oc S{y). 
This yields lmi;^'"(0) cx T. On the other hand, at 
higher temperatures (O.lfis < T < il^ in Fig. 9, re- 
ferred to herein as region II), the strength of lm^„m(0) 
is weak and the temperature dependence comes from 
~ 77B(f2^)Im^„'m' (rj'g), where Q'^ is the renormalized 
Einstein energy. Although Q'^ is weakly temperature 
dependent, for simplicity, we regard fl'^ as a temper- 
ature independent parameter. These observations yield 
ImE^'"(0) oc nein'^) ^ e^pi-fl'^/T), for T < n'^. 

For the phonon Green function, eq. (A-4) is reduced, 
at low temperature and frequency with z <C T to 



X ImAjnn'iy + z)linArn'„{y), 



where we have used 



(41) 
1 for 



d^e-^/^mAuix) 

T < Qe- To simplify eq. (41) further, we omit the in- 
dices 71, m, n' , and to'. Taking into account the fact 
that the dominant contribution comes from y = in 
the integrand, and -lmA(y) ~ {^/2)/{y'^ + (7/2)^) with 
7 = — 2ImS^(0) (this is of course too simplified an ap- 
proximation), we can estimate —lmD{z)/z as 



lmD{z) 



1 



oc 



7 



.0 ^ 
"7^7' 



oc 



T z2 + 72 

T-2 
r-iexp(f]^/r) 

These estimations reproduce the numerical results shown 
in Fig. 9. The expression eq. (42) is expected to be a rele- 



in I, 
in II. 



(42) 
(43) 
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Fig. 10. Temperature dependence of ReDio„{z) with different 
frequencies with A = 0.35Qe = 0.0035. The parameters are as 
in Fig. 5: g = 0.01, A^cut = 6 and v = 0.1. 

vant origin for the Debye-type formula in eq. (1). Indeed, 
the form of eq. (42) is as same as an expression in the 
case in which there are degenerate modes with the life- 
time T = 'y~^ in the ground state.'^^ This is manifest by 
the 1/T in front of eq. (42). 

3.3 Frequency dependence on sound velocity 

In order to estimate the frequency dependence on the 
sound velocity correction with the use of the Holstein- 
Anderson model, we note that the high-energy spectra 
ImD(z) does not contribute to the frequency depen- 
dence in ReD{z) for z ^ fis through the Kramers- 
Kronig transformation. The high-frequency part only 
gives weak monotonic temperature dependence. Thus, 
we define ReZ?iow(2;) by introducing the cutoff parameter 
A < as: 

/ s ^ dxlmD(x) , , 

ReAow(^) = P / —. (44) 

J_A TT X~ Z 

From the discussions in §2, we can estimate the effect 
of the optical phonon on the sound velocity. We approxi- 
mately replace the second term in eq. (14) by the phonon 
Green function obtained in the previous single-site prob- 
lem: 

u\zA)^ul{q)^\^^D{z). (45) 

Although the present results cannot be directly com- 
pared to the experimental data in the case in which the 
parity of the optical mode is —1, it is expected that their 
qualitative behaviors do not differ greatly. 

In Fig. 10, we show the temperature dependence of 
ReDiow(2;) with different frequencies with A = 0.0035 — 
0.35 X rig and g — ^e- Although the choice of A includes 
a certain arbitrariness, the essential feature of ReZ^iow 
does not change if we take A to be smaller than ~ O.SfJe. 
In the case of harmonic coupling, Re£'iow(^) is directly 
related to the sound velocity and the EC (we cannot 
discuss the anisotropy in EC's, see §4). For this case, 



the frequency-dependent term C(z) in the EC's is simply 
given by 

C(z) cxReAow(z). (46) 

Clear hardenings, the positions of which shift to higher 
temperatures as the frequency increases, can be observed 
in Fig. 10. This feature is consistent with the phenomeno- 
logical expression (1). 

In addition to ReDiow , there are contributions from the 
high-energy spectra: Re£'high(2) = Re£'(z)— Re£'iow(2) 
and the anharmonicity of acoustic phonons, both of 
which are essentially frequency independent. The elastic 
constant observed in the experiments is basically given 
by these two contributions unless we take into account 
electric contributions. The former makes the elastic con- 
stant soften at low temperatures and low frequency. On 
the other hand, the latter makes the elastic constant 
harden. We show the temperature, g (coupling constant), 
and frequency dependence of ReDiow and ReDhigh in 
Figs. ll(a)-ll(d). The magnitude of ReZ^iow increases 
as g increases compared to the ReUhigh- Although, in 
Fig. 11(a), there exist shoulders as a function of temper- 
atures in Re£'iow(-z), as shown in the inset of Fig. 11(a), 
their magnitudes are quite small relative to Rel^high- On 
the other hand, these features disappear at larger g, as 
shown in Fig. 11(c). The shoulder in the parameters in 
Fig. 11(c) can be seen at the higher frequency (not shown 
in Fig. 11(c)). In all of the data in Figs. 11(a)- 11(c), soft- 
ening behaviors are observed at ReZ?iow(-z) with small z. 
These features originate from the development of the soft 
mode being too strong in our calculation. More reliable 
methods are needed for the low-temperature region in 
order to discuss this point further. 

In order to observe the temperature dependence of 
the shoulder in ReD(z), it is useful to differentiate 
ReDiow(z) by T . In Fig. 12, we plot the temperatures T*, 
where \dBjs^D\a^(z) j dT\ takes a maximum value. Roughly 
speaking, T* corresponds to the midpoint of the shoul- 
der. In Fig. 12, we show the results for three different 
parameter sets, keeping the renormalized Einstein fre- 
quency Sl'g the same for each set. As shown in Fig. 12, 
the frequency dependence of fi^/T* is well described by 
the condition zt = 1 with r = Toexp(a"/T). Here, as 
before, a" ~ 0.717^ is a parameter on the order of fi^. 
Interestingly, the estimated values of are approxi- 
mately a" X 1/100. This is in good qualitative agreement 
with the experimental values Tq^ = O.OOGa" — O.la" 
for filled skutterudites and clathrate compounds. We can 
find that Tq^ becomes small as g decreases. This arises 
from the fact that the soft mode amplitude becomes 
smaller and its frequency region moves to the lower side 
as the electron-phonon interaction g becomes small. This 
feature is sensitive to the activation energy a". A small 
change in a" can affect the magnitude of tq. This is the 
reason why we show the data with the same il^ values 
in Fig. 12. 

4. Discussions 

Thus far, we have restricted our study to the examina- 
tion of a system with one component phonon. Here, we 
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discuss qualitatively the anisotropy in the UD observed 
in the above experiments. 

In La3Pd2oGe5, anomalous temperature dependence in 
Raman spectra are reported in the "T2(," mode of La.^° 
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Fig. 12. Frequency dependence of T* for three different g values, 
maintaining the renormalized Einstein frequency £7^ constant. 
f2'g is determined by the peak position of lmD{z) at the high- 
energy region in the zero temperature limit. We plot the dashed 
lines, which are determined by the condition ZToexp(a"/T) = 
1. The values of ro and a" are listed in the panel. The other 
parameters are as in Fig. 5. In particular, the data for g/O^ = 
1.14 is the same as in Fig. 5. 



For this material, our results suggest that this mode 
strongly interacts with conduction electrons, and the ob- 
served UD in C44 can be naturally interpreted according 
to the discussion in §2.1. For the optical phonon of the 
rare-earth dominated mode (Ti„) in filled-skutterudites, 
there are three degrees of freedom, namely x, y, and z. 
In order to discuss the anisotropy in the UD, it is desired 
to take into account these aspects of the system. For ex- 
ample, we must estimate the wave number dependence 
in the electron-phonon coupling g, which has a crucial 
role for the anisotropy in the UD. The results obtained 
in the present paper can be thought of as simplified but 
essential for describing the frequency dependence in the 
EC. Experimentally, there is an inconsistent result re- 
garding the anisotropy of the UD in the different samples 
of PrOs4Sbi2.^^ This discrepancy should be investigated 
by further experimental efforts. 

At this stage, the question as to whether the rare-earth 
Tiu mode is the relevant mode for the realization of the 
UD in filled skutterudites arises. It might be interest- 
ing to determine whether the temperature dependence 
of the Eg modes in filled skutterudites is unusual, e.g., 
whether softening occurs. It is important to carry out 
high-resolution Raman scattering experiments to clarify 
this point. Recently, Ogita et ai, observed the second- 
order Raman spectra in filled skutterudites (especially 
Sb compounds) as an anomalous property of the rare- 
earth mode.^^ In the case in which the relevant mode is 
actually the rare-earth Ti„ mode in filled skutterudites, 
the mechanism of the anisotropy in the UD might be 
more complicated, as mentioned above. The Ti„ opti- 
cal phonon should be investigated as a relevant optical 
phonon for the UD as a further theoretical study. Ex- 
periments in unfilled skutterudites with electric bands 
similar to, e.g., LaOs4Sbi2 are also desired in order to 
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identify which mode is important for the UD. Applying 
this theory to optical phonons with even parity and com- 
paring the energies of the optical phonons estimated by 
Raman scattering and the activation energy obtained by 
the ultrasound experiments for each materials will pro- 
vide various interesting properties. 

In filled skutterudites, anomalous phonon contribu- 
tions are observed in various quantities, especially in 
ROs4Sbi2 (R=Pr, La, Sm, etc.). One reason why the 
effect is prominent in ROs4Sbi2 may be that the size of 
the Sbi2 cage is the largest among existing filled skut- 
terudites, including the rare-earth ion R. Another rea- 
son may be that the conduction electrons near the Fermi 
level in Os-compounds consist of a molecular orbital of 
the cage (Ai„) and the d-electron in Os. Unlike the Ru- 
compound, in which the levels of 4d electrons are deep 
compared to that of 5d or 3d electrons, so that their 
components near the Fermi level are negligible, the con- 
tributions from the d-electrons are much larger. It is thus 
expected that d-electrons play an important role in Os 
compounds. The latter point, however, requires more so- 
phisticated analysis, because Iwasa et al., recently ob- 
served the quasi-elastic peak even in Ru compounds. 
The observed width of the quasi-elastic peak is too large 
for UD to be realized in the MHz frequency range. In or- 
der to clarify the reason why the UD has been observed 
in Os-compounds but not in Ru-compounds, it may be 
necessary to carry out experiments in a wider frequency 
range. Note that if we assume the dominant electron- 
phonon coupling between the Aiu molecular orbital and 
the T\u optical phonon is local and linear in the displace- 
ments, it is impossible for the A\u molecular orbital to 
interact with the Ti„ optical phonon due to symmetry. 

Recently, nuclear magnetic resonance (NMR) experi- 
ments on KOs206^^ and LaOs4Sbi2^^ were carried out. 
The results indicate that the relaxation is due to the 
quadrupolar coupling between the nuclei and the ionic 
motion (through a direct process) at K and La sites. In 
this case, l/(TiT) is given by a phonon spectral weight:^^ 
l/(TiT) (X — ImD(z)/z. The results are interpreted by a 
strongly damped oscillator with an activation-type relax- 
ation time. Although the energy scale of the ultrasonic 
and NMR experiments is of the same order as the magni- 
tude (~ MHz) , the results show no frequency dependence 
in l/(TiT). This might be due to the difference between 
the q = and g-integrated spectrmn, or the difference 
of the phonon mode in the ultrasonic and NMR exper- 
iments, as mentioned above (the ultrasonic experiment 
for KOS2O6 has not yet been performed). Based on Fig. 
6, it is expected that l/(TiT) has frequency dependence. 
In the case in which the relaxation time is indeed an 
activation-type, i.e., oc cxp{E^/T), the small change in 
z does not indicate the same order change occurs in T, 
due to the exponential dependence. The precise origin of 
this problem is not yet understood and further experi- 
mental and theoretical studies are required. 

For the theoretical aspects, the approximated method 
used in the present paper becomes worse in the low- 
temperature regime, because —lim.z^o^T^D{z)/z di- 
verges as T~^. This might cause the unphysical develop- 
ment of the soft mode. Intuitively, the divergence must 



be stopped at a certain temperature. However, it is im- 
portant for this method to be able to capture the tem- 
perature dependence of the —lmD{z) for z <C at the 
intermediate temperature regions, where UD is observed. 

In our treatment of pseudo-fermion representations, it 
is straightforward to use higher-order interactions e.g., 
X"^, X^, ■■■ as interaction parts of the Hamiltonian. 
This is because the same argument can be applicable 
by simply replacing the matrix elements X with X". 
The model including these higher-order electron-plionon 
couplings becomes important for ions with very large 
displacements, such as the rare-earth ion in filled skut- 
terudites. The variations of hybridization between the f- 
electron and conduction electron by these ionic motions 
should be taken into account in a realistic model for the 
further clarification of electron-phonon systems. For the 
anharmonicity of the potential for the local phonon, it 
is straightforward to include higher-order terms such as 
X^. Although the same line of discussion can be intro- 
duced in the present paper from an off-center potential, 
the presence of the off-center potential does not play a 
fundamental role with respect to our results. 

Finally, we mention the relationship between the 
present theory and our previous studies. In the 
present paper, we have discussed primarily the proper- 
ties of a local phonon that interacts with conduction elec- 
trons. In our previous studies, however, we discussed the 
properties of electrons interacting with off-center config- 
urations of an ion, which is an effective theory for the 
strong coupling limit of an electron-phonon system. 
Further theoretical studies are needed in order to unify 
these aspects. 

5. Summary 

In conclusion, we have investigated the acoustic and 
optical phonon spectra in a strongly coupled electron- 
phonon system at finite temperatures. We have ap- 
plied the self-consistent ladder approximation, which was 
used successfully in the Coqblin-Schrieffer model with 
crystalline-electric-field states of Ce impurities, to the 
electron-phonon system, and have determined that the 
low-energy peak in the optical phonon spectral function 
develops at low-tcmpcrat\ires and clauses the ultrasonic 
dispersion observed in fiUed-skutterudites and clathrate 
compounds. The temperature dependence of this peak 
shows activation type behaviors at the intermediate tem- 
perature regions, leading to the Debye-type formula for 
the ultrasonic dispersion with activation-type relaxation 
time. These results qualitatively explain the experimen- 
tal results of the ultrasonic region observed in the filled 
skutterudites and related compounds. 
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Appendix: Self-consistent Equations for the 
Self-Energies and Phonon Green 
Function in a Real Frequency 

In the Appendix, we list the self-consistent equation 
for the imaginary part of the retarded self-energies and 
the form of the retarded phonon Green function. 

Self-Energies: 



dy 



nF{y)lTxiG{-y) 



[ J2 Xnm'lmxdA {y + z)Xn'm (A-1) 



and 



ImSG(z) = (1 



n 



E 



dx - n-l 



with 



dy 



(A2) 



-nF{y)lmG{y)lmAn^{z + 2/)(A-3) 



where xga in eq. (A- 2) is obtained using eq. (30). 
In eqs. (A-1). (A-2). and (A-3), np represents the 
Fermi distribution function. The function with an 
overbar is defined as f{x) = f{x + X), taking A — > c» 
subject to the local constraint '^n'^" ^ 1. In cq. 
(A-2), we have divided the result, which is obtained 
by the analytic continuation of eq. (28), by X^„(aJiCin)-^^ 



Phonon Green Function: 
lmD{z) = -(1-e-^ 



X ImAmn'iy + z)lTnAm>„{y)X„mXn'm' 

X [-E/ f e-/^W„(x)]"\ (A-4) 

In eq. (A-4), we have used the same prescription as that 
used for ImSg. 
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